Dark radiation as a signature of dark energy 
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We propose a simple dark energy model with the following properties: the model predicts a 
late-time dark radiation component that is not ruled out by current observational data, but which 
produces a distinctive time-dependent equation of state w(z) for z < 3. The dark energy field can 
be coupled strongly enough to Standard Model particles to be detected in colliders, and the model 
requires only modest additional particle content and little or no fine-tuning other than a new energy 
scale of order milli-electron volts. 
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I. INTRODUCTION 

Considerable evidence [l|, HJ has accumulated suggest- 
ing that approximately 70% of the energy density in 
the Universe comes in the form of an exotic, negative- 
pressure component, called dark energy. (For a recent 
review, see Q.) 

The equation of state (EoS) parameter is defined as 
the ratio of the dark energy pressure to its density: 



w 



Pde/pde • 



(1) 



Observations constrain w to be very close to — 1. For 
instance, if w is assumed to be constant, then necessarily 
-1.1 < w < - 0.9 |, |. If w = -1, the dark en- 
ergy density remains constant even though the Universe 
is expanding. The simplest way of producing a w = — 1 
component is through a cosmological constant, or vac- 
uum energy density. However, as is well known, the en- 
ergy density needed to explain the observed acceleration, 
A 4 = (lO _3 eV) , is considerably smaller than the value 

of (10 19 GeV) 4 (Planck density) predicted from quantum 
field theory This 124-orders-of-magnitude discrepancy is 
called the cosmological constant problem. 

The fact that the observed vacuum energy also hap- 
pens to be just a few times greater than the present mat- 
ter density has led to speculations that it might in fact 
be evolving with time - only now reaching a value com- 
parable to the matter density. Such a time-varying vac- 
uum energy is sometimes referred to as quintessence. The 
simplest way of achieving a time- varying vacuum energy 
is through the use of s pat ially homogeneous canonical 
scalar fields @, 0, H, [9l. TToj. In these models, the field 
typically rolls down a very shallow potential, eventually 
coming to rest when it can find a local minimum. 

Quintessence models typically have fine-tuning prob- 
lems. For example, since the quintessence redshifts 



more slowly than ordinary matter or radiation, the cur- 
rent quintessence dominance can only be explained by 
fine-tuning the initial conditions. This problem can be 
avoided in the class of so-called "tracker" models, in 
which the evolution of the quintessence field is insensitive 
to the initial conditions. For generic quintessence mod- 
els, the flatness of the potential makes any excitations 
of the field almost massless ~ 10~ 33 eV. To provide the 
necessary vacuum energy density, the present value of the 
potential energy should be on the order of A 4 (although 
there is really no rigorous physical reason to expect this). 
The field value 0o today should therefore be on the or- 
der of the Planck mass, i.e. (j> - lO^GeV. 1 In [l7[ 
it was shown that couplings between quintessence and 
ordinary matter, even if Planck-suppressed, can lead to 
long range forces and time-dependence in the constants of 
Nature, both of which are tightly constrained. Reference 
[Tc| showed that even Planck-suppressed thermal interac- 
tions between matter and quintessence can significantly 
alter the evolution of the latter, leading to a problematic 
equation of state. 

An alternative to a slowly rolling field is a scenario 
where the field is stuck in a false vacuum minimum. 
In this case, the observed cosmological constant is at- 
tributed to the energy difference between the false and 
true vacua, which could either arise from higher-order 
or non-perturbative effects [2l[ . Other proposals 
for the origin of the false vacuum energy are the con- 
fining scale of a hidden SU(2) sector |22j, Planck-scale 
suppressed mediation into a hidden sector of electroweak 
TeV-scale supersymmetry breaking j !23l . |24| , or the vac- 
uum energy of a hidden sector which is stuck in a state 
of equilibrium between phases [25[ . 

In many of these quintessence models, the field(s) re- 
sponsible for the acceleration have to be almost com- 
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1 An alternative class of models which relies on non-linear oscilla- 
tions of the quintessence field and does not require an extremely 
flat potential was proposed in 111! I l'J and discussed further in 
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pletely decoupled from the rest of the Universe. 2 This 
is disappointing, since it suggests that direct detection 
of quintessence through its interactions with Standard 
Model particles will be extremely challenging, perhaps 
impossible. 

In this paper we present a quintessence scenario in 
which the dark energy field can be coupled strongly 
enough to Standard Model particles to be detected in col- 
liders, and which allows for a significant time variation 
in the equation of state. This time-varying w = w(z) 
has a characteristic form which depends on only a single 
parameter, and can thus be excluded by cosmological ob- 
servations in the near future. Our model only requires a 
singlet scalar field (or, alternatively, a small gauge sector 
like SU (3) Yang-Mills theory; other possible realizations 
are also briefly outlined at the end of the paper) and a 
new energy scale on the order of milli-electron volts. 

II. MODEL 

Consider a singlet scalar field dark energy with La- 
grangian 

C= l -(d,<t>f~V{<t>). (2) 

We allow this field to be strongly coupled to Standard 
Model particles. The finite temperature effective poten- 
tial, which includes interactions of this field with virtual 
particles and the heat bath, can be taken to be similar to 
the Higgs potential in the electroweak phase transition 
(see, e.g., [27| for a review): 

V (4>,T) = A + D (T 2 - T 2 2 ) (f> 2 - ETcf) 3 + ^A^ 4 . (3) 

D, E, A and A are constants. A can be adjusted to 
give the correct value of the observed dark energy den- 
sity when T = 0. T 2 is defined as the temperature where 
V" (<j> = 0) = 0. We choose T 2 = A, i.e. roughly T 2 ~ 
11.6 K, and assume that it represents a new energy scale 
in particle physics. At high temperatures, T 3> T%, <p — 
is the only minimum of the potential. As the Universe 
cools down, an inflection point appears in the potential at 
temperature T* = T 2 /y/l - 9E 2 /8XD. At lower temper- 
atures, this splits into a barrier and a second minimum. 
The critical temperature Ti = T 2 /y/l - E 2 /\D corre- 
sponds to the point where the second minimum is equal 
in (free) energy to the <p — minimum. At temperatures 
T < T\, the second minimum has lower free energy than 
the one at <f> = 0. The evolution of the potential well 
with temperature is shown in Fig. ((T|). 



2 For an entirely different type of dark energy model, which can 
have a particle physics signature, see |2(I . Here the acceleration 
is provided by mass-varying neutrinos (MaVaNs) which act as a 
negative-pressure fluid. 
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FIG. 1: An example for the evolution of the finite temperature 
effective potential V {4>, T) of the dark energy field, eqn. 
as the temperature T decreases through the first order phase 
transition region ~ Ti. 

Let us now consider the dynamics of the dark energy 
field. At temperatures T > T\ (which corresponds to 
roughly z > 3) the dark energy field remains trapped at 
the <f> = minimum, providing a constant energy density, 
which we assume to be slightly higher than A 4 . As the 
temperature approaches T± and below, a first order phase 
transition is triggered as the field tunnels into the true 
vacuum A. The physics of the phase transition is almost 
identical to that of the Higgs sector in models of elec- 
troweak baryogenesis. This transition releases energy in 
relativistic modes (i.e., scalar particles of the cf> field), and 
brings the vacuum energy to A. Because the correlation 
length of the transition is microscopic, and the relativis- 
tic modes couple weakly to ordinary matter (i.e., more 
weakly than photons, perhaps similar to neutrinos), such 
a transition is only loosely constrained by observation. 
The positive pressure of the radiation, which eventually 
redshifts away, causes the effective EoS of the dark energy 
to vary in (redshift) time z. 

We note that the only important feature of the model 
described above is that it has a weakly first order phase 
transition at a temperature of order A, which is natural if 
one assumes the dynamics of (f> to be entirely determined 
by that energy scale and dimensionless couplings of order 
one. It is an interesting coincidence that this occurs at 
a redshift of z ~ 3 if the temperature of the dark energy 
field is similar to that of the Standard Model particles. 
This need not be the case, but it seems a reasonable 
assumption, especially if there are non-negligible inter- 
actions between <f> and ordinary particles, which would 
enforce thermal equilibrium at sufficiently high tempera- 
tures. When the transition happens at z ~ 3 the result- 
ing radiation component leads to significant and charac- 
teristic variation in w(z). The form of w(z) is determined 
by a single parameter - the energy fraction in relativistic 
dark radiation modes just after the phase transition. In 
some cases, such as the gauge models discussed below, 
even this fraction is calculable from the phase diagram. 

Any sector which produces a weakly first order tran- 



3 



sition at a temperature of order A would also suffice. 
For example, pure SU(N) gauge theories with N > 2 
have first order deconfinement phase transitions [28| and 
exhibit effective potentials like those in Fig. (Q]), with 
(j> an order parameter for confinement, for example the 
Polyakov loop. Here, the latent heat and fraction of en- 
ergy in relativistic modes is calculable via lattice simula- 
tion. 

We stress that the models discussed do not in any way 
explain the existence of the energy scale A, or why it 
determines the vacuum energy density today. In par- 
ticular, why should the vacuum energies from all the 
other degrees of freedom cancel out, leaving the dark en- 
ergy field to determine the cosmological constant? One 
way of explaining this would be to assume that some- 
where in the configuration space, outside the region de- 
picted in Fig. {l}, t ne potential reaches a global mini- 
mum V (</>*) = 0, where the total vacuum energy (in- 
cluding zero point energies and radiative corrections from 
all fields) is exactly zero. That is, some currently un- 
known mechanism (Euclidean wormholes, quantum grav- 
ity, ...) conspires to make the total vacuum energy zero 
at (f> = 0*, implying that the deviation of V (</>) from zero 
is the only vacuum energy. 

In any case, if one assumes that new physics at the 
energy scale A determines the observed cosmological 
constant, it is easy to obtain a predictable redshift- 
dependent w — w(z) together with interesting particle 
physics signatures - no fine tuning of parameters is re- 
quired. 



III. OBSERVATIONAL CONSEQUENCES 
A. Astrophysics 

As discussed in the previous section, our model pro- 
duces a certain amount of dark radiation at redshift 
z ~ 3. This radiation affects the Hubble expansion rate 
H as well as the effective equation of state of the dark 
energy. 

Let / denote the fraction of the dark energy that is to- 
day in the form of relativistic modes, and let us assume 
that the phase transition occurred at redshift zpt = 3. 
The Hubble expansion rate H (z) after the phase transi- 
tion z < zpt can be written as 



H 2 (z) = H 2 Q m0 (1 + zf + Q r0 (1 + zf + 



+fn^ (i + zf + n 4 



(4) 



where Q m o, ^ro and f^o denote the present-day values 
of the density parameters of matter, radiation and dark 
energy. Note that at sufficiently low temperatures the 
non-zero mass of the dark radiation (coming, e.g., from 
the curvature at the lower minimum in Fig. (TT|)) will be 
non-negligible and its energy density will then redshift as 
(1 + z) 3 . 
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FIG. 2: Likelihood contour for the parameters / and Q m o- 
The yellow (light) region is excluded at the 2a level and the 
orange (darker) region is excluded at the la level. The red 
(darkest) region is not excluded at either confidence level. 



Our model superficially resembles other models with 
a dark radiation component, such as models with extra 
relativistic degrees of freedom or the Randall-Sundrum 
model with dark radiation. The difference, of course, is 
that in our model the dark radiation arises very late, and 
so is not subject to the well-known limits from Big Bang 
nucleosynthesis or the cosmic microwave background. It 
was noted by Zentner and Walker [2^] that if one consid- 
ers only late-time constraints on extra relativistic degrees 
of freedom from SNIa data, the limits are surprisingly 
weak. Our results, which we describe now, agree with 
this conclusion, even with the addition of more recent 
SNIa data. 

In Fig. ([2|) we construct a likelihood plot for the param- 
eters fi m o and /. We choose fi^o = 0.7 and marginalize 
over the present value of the Hubble parameter Hq us- 
ing the recent Type la Supernovae standard candle data 
(ESSENCE+SNLS+HST) from [5]. Clearly, the SNIa 
data do not rule out a sizable fraction of the dark energy 
today being in relativistic modes. 

It is easy to derive an analytic expression for w(z) in 
this model. Taking pdr, Pdr to be the dark radiation 
pressure and density, and p^, p^ to be the scalar field 
pressure and density, we have w = (pdr+P<p)/ (pdr+P4>)- 
But p D R — pdr/ ^ and p^ = —p^, leading to 



w(z) 



(l/3)/(l + z) 4 -l 
/(I + zY + 1 



(5) 



As an example of the possible strong late-time variation 
of w(z) predicted by this model, in Fig. ^ we plot w 
vs. z for / = 0.01; e.g., the parameter choice / = 0.01, 
!l m (] = 0.28 is conservative, it is not excluded at la by 
the SNIa data. The pressure of the relativistic component 
increases the effective EoS of the dark energy component 
with increasing redshift. Note that, in our model, the 
shape of the w vs. z curve is fixed once / is fixed. 

Another diagnostic for dark energy models is the evolu- 
tion in the w — w' phase plane (30| . where w' = a(dw/da) 
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FIG. 3: w vs. 2 for the choice / = 0.01, which along with, 
e.g., a conservative Q m o = 0.28 is not excluded at the la level 
by the SNIa data (see Fig. ©). 

and a = 1/(1 + 2) is the scale factor. From (O, we find 
that 
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FIG. 4: Likelihood contour for the parameters fi m o and z*, 
the redshift at which the Universe has exited the false vacuum 
and entered the true vacuum, releasing the energy of the cos- 
mological constant into relativistic modes. The yellow (light) 
region is excluded at the 2a level, and the orange (darker) 
region is excluded at the la level. The red (darkest) region is 
not excluded at either confidence level. 



w 



(l + w)(3w - 1) 



(6) 



Note that the relationship between w' and w is indepen- 
dent of /. This means that all of these models evolve 
along the same evolutionary track in the w — w' plane; 
the value of / simply determines where the model sits on 
this evolutionary path at the present time. 

In the terminology of [3fJ, these are "freezing" 
quintessence models, since w decreases with time to 
w — > — 1. However, eqn. © predicts behavior that is 
distinct from standard freezing quintessence models. For 
them, [30j suggested the bound w' > 3w(l + w), whereas 
our model always has w' < 3w(l + w). In this respect, it 
more cl osely resembles the barotropic models discussed 
in refs. [3l|, [32]. This result arises from the fact that 
we have a two-component dark energy model. In terms 
of the evolution of the equation of state, our model re- 
sembles the barotropic "wet dark fluid" model proposed 
in [33|, HH , with the important difference that the dark 
radiation in our model appears only at late times. 

As a point related to our analysis, we consider the 
possibility that our Universe has exited the false vacuum 
in recent times, i.e. all of the dark energy has recently 
been dumped into relativistic modes, which will eventu- 
ally redshift away. This would be the case if the dark 
energy field went through a first-order phase transition 
of the type considered above, but into the true V = 
vacuum and not into another meta-stable vacuum. This 
scenario also arises in the "accelerescence" model consid- 
ered in [24j. Let z* be the redshift-time of this phase 
transition, when the vacuum energy is instantaneously 
(relative to cosmological timescales) converted into radi- 
ation. The Hubble parameter is therefore given by 



h 2 (z) = Hi n m0 (1 + z y + n r0 (1 + zf + 



-l-fi^o (1 + zS {z* - z)) 1 



where S (x) is the Heaviside step function. Using the 
SNIa data, Fig. ((4| is a likelihood plot for the parameters 
Q m ,o and z*. We find that z* is tightly constrained by 
the data, with the maximum z* allowed being ~ 0.1 at 
2a. Thus, if the Universe has already exited the vacuum 
energy epoch, it did so very recently. 

B. Particle physics 

An interesting feature of our scenario is that the dark 
energy field can be coupled relatively strongly to Stan- 
dard Model particles. This makes it possible, in principle, 
for this kind of dark energy to be detected in colliders. 3 



(7) 



3 The proposed dark energy field <f> (dark radiation) has very 
small mass ~ meV, and might be produced by thermal reac- 
tions in stars. The energy loss argument for globular-cluster 
stars or red giants sets some strict limits on its coupling to the 
Standard Model, similar to constraints on the axion decay con- 
stant |35(| . E.g., an Oil) Yukawa coupling of a scalar field <j> to 
quarks q induces, in one-loop, an effective dimension-5 coupling 
to photons a<j> {F^) 2 /4irm q ; this coupling would cause glob- 
ular clusters to lose energy (into <f> modes) more quickly than 
is actually observed, unless suppressed by a quark mass scale 
m q > 10 7 GeV, thus disallowing such a coupling to any Standard 
Model fermions. On the other hand, if the dark energy field obeys 
a Z2 symmetry <fi — * — tfi, or if tj> is the glueball field of some addi- 
tional SU(N) gauge theory (with interpolating dimension-4 op- 
erator G M „G M ") coupled to the Standard Model via messengers 
m q , the induced effective interaction with photons has higher 
dimension; dimension-6 interactions a<f> 2 (F^) 2 /Airin 2 can al- 
ready avoid the energy-loss constraints for helium-burning stars 
(T core ~ 10 s K) if m,j<; 20 GeV, thereby allowing coupling of 
the dark energy field to weak-scale Standard Model particles. 
These astrophysical constraints do not significantly hinder de- 
tection of our proposed dark energy field <f> at particle colliders, 
which provide energies 2> T corc and produce weak-scale particles 
abundantly. 
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FIG. 5: Potential energy surface for a gauge theory, where 
<j> is an order parameter for confinement (Polyakov loop) and 
<f> a colored scalar field. For N > 2, SU(N) models will 
have a first order confinement-deconfinement transition as the 
temperature is lowered. However, at zero temperature the 
deconfined phase ($ = (f> — 0) is not necessarily metastable. 



The simplest model we considered, comprised of a sin- 
glet scalar (f>, has some challeng direct coupling 
between cj> and the Higgs boson operator WH cannot be 
excluded. This would lead to significant radiative correc- 
tions to the 4> potential parameters, making the model 
somewhat unnatural. However, if this fine tuning is ig- 
nored, the <p-H coupling would provide for direct pro- 
duction of (j> particles at colliders. 

Our alternative model uses a pure SU (N) gauge theory 
sector (N > 2) with strong coupling scale A ~ A. This 
model requires no fine tuning and the fraction of energy 
in relativistic modes after the phase transition can in 
principle be calculated from simulations of the SU(N) 
theory. Glueballs of this sector would be light excitations 
with mass of order A; the phase transition temperature 
would be at least a few times the glueball mass. The 
glueballs could couple to Standard Model particles via 
higher dimension operators such as 

(*fiv ftm ) (8) 

where G is the SU(N) field strength and O sm a (Lorentz 
scalar) Standard Model operator such as H^H, qq, etc. 

If we wish to ensure that there exists a point in the 
configuration space where the vacuum energy vanishes 
exactly, one must add some extra degrees of freedom. For 
example, a colored scalar $, whose potential V ($) has 
positive second derivative at $ = and a global minimum 
at some non-zero value, would suffice (see Fig. ©). Note, 
this likely requires a non-renormalizable potential (i.e., 
with <£> 6 term). 

An exactly vanishing potential energy at some point 
in the configuration space is a generic feature of many 
theories with global supersymmetry [3(| - the vacuum 
energy is zero precisely at the supersymmetric points. 
While this fact does not explain away the cosmological 
constant term in the Einstein-Hilbert action, it may have 



something to do with the existence of an absolute mini- 
mum with small or vanishing energy density. In the su- 
persymmetric framework, a presently non-zero and posi- 
tive vacuum energy can be explained by the fact that the 
Universe is currently sitting at a meta-stable vacuum of 
the field theory, and its difference to a supersymmetric 
vacuum gives the present positive vacuum energy den- 
sity A 4 . Examples of supersymmetric theories with such 
meta-stable vacua can readily be given either as simple 
Wess-Zumino models (e.g. [37j ) or in terms of supersym- 
metric gauge theories which provide an ultraviolet frame- 
work for (O'Raifeartaigh-like) meta-stable supersymme- 
try breaking [38] . Furthermore, the dynamics for the 
Universe to initially be stuck in a meta-stable vacuum 
in the course of its cooling with only subsequent transi- 
tion to the absolute supersymmetric minimum has been 
confirmed [39[ and the lifetime of such meta-stable vacua 
has been considered. 

The phenomenologically plausible scale for electroweak 
supersymmetry breaking of 1 TeV 3> A is much too large 
to directly account for the observed vacuum energy in the 
way just outlined. Nevertheless, supersymmetry might 
be invoked to provide for an absolute zero of the energy, 
at least in the sector containing the dark energy dynam- 
ics itself. By small modifications to the toy models de- 
scribed in the previous paragraph (changing parameters, 
or adding one or two new chiral superfields) it is further- 
more possible to build dark energy sectors, containing 
just a small number of chiral superfields (possibly arising 
as effective fields [38j]), with two slightly non-degenerate 
meta-stable vacua of energy ~ A along with a supersym- 
metric minimum, thereby giving a natural explanation 
for an absolute zero energy and also exhibiting the inter- 
esting dynamics of dark energy and recent dark radiation 
of our models @. 

The dark energy sector is likely to feel electroweak su- 
persymmetry breaking, at least through gravitational ef- 
fects [24j], and is therefore not expected to be perfectly 
supersymmetric. If this mediation happens only through 
gravitational interactions, the terms induced in the dark 
energy sector are naturally of the correct order of magni- 
tude (1 TeV) 2 /Mpi ~ 10 -3 eV ~ A to provide for energy 
differences in the dark sector of the size of the observed 
cosmological constant. Furthermore, a slight modifica- 
tion of the first model of [24J, e.g. addition of an m$ 2 
term to the superpotential of the dark sector, generically 
yields three non-degenerate meta-stable vacua of energy 
~ m ~ A, again yielding our scenario of dark radiation 
along with dark energy. If it is implemented in Nature, 
there could be dark radiation according to |(3J) as well as 
dark energy of order A present. And unlike in the model 
in [2^|, this scenario would cosmologically be detectable 
not solely in the very far future (billions of years from 
now), but could be confirmed, rejected or constrained al- 
ready in the foreseeable future through comparison of the 
more precisely measured past expansion rate of the Uni- 
verse (for 0.1 < z < 3) to our predictions for the equation 
of state w(z) as Fig. 
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IV. CONCLUSIONS 

We have discussed a class of dark energy models which 
have interesting cosmological as well as collider signa- 
tures. In these models, a first-order phase transition at 
redshift z ~ 3 releases energy in relativistic modes (dark 
radiation) leading to a characteristic time-dependence in 
the effective dark energy equation of state. We have 
shown that such models are consistent with SNIa data, 
and are relatively easy to construct as extensions to the 
Standard Model. 

As an interesting and important side issue, we consid- 
ered the possibility that the Universe might have recently 
(at redshift z = Z*) exited the false vacuum phase and 



entered the true vacuum, converting all of the dark en- 
ergy into relativistic modes. We show that the SNIa data 
places tight constraints on z„, restricting it to z* ~ 0.1 
or less at the 2a confidence level. 
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